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Boundary layer receptivity and control 


By D. C. Hill 


1. Motivation and objectives 

Receptivity processes initiate natural instabilities in a boundary layer. The in- 
stabilities grow and eventually break down to turbulence. Consequently, receptivity 
questions are a critical element of the analysis of the transition process. Success 
in modeling the physics of receptivity processes thus has a direct bearing on tech- 
nological issues of drag reduction. The means by which transitional flows can be 
controlled is also a major concern: questions of control are tied inevitably to those 
of receptivity. 

Adjoint systems provide a highly effective mathematical method for approaching 
many of the questions associated with both receptivity and control. The reader 
is referred to Hill (1993) for a detailed description of their use in the receptivity 
context. The long term objective of this project is to develop adjoint methods to 
handle increasingly complex receptivity questions, and to find systematic procedures 
for deducing effective control strategies. 

The most elementary receptivity problem is that in which a parallel boundary 
layer is forced by time-harmonic sources of various types. The characteristics of 
the response to such forcing form the building blocks for more complex receptivity 
mechanisms. The first objective of this year’s research effort was to investigate how 
a parallel Blasius boundary layer responds to general direct forcing. 

Acoustic disturbances in the freestream can be scattered by flow non- uniformities 
to produce Tollmien-Schlichting waves. For example, scattering by surface rough- 
ness is known to provide an efficient receptivity path. This problem has been 
investigated previously in a number of different ways. The present effort is directed 
towards finding a solution by a simple adjoint analysis, because adjoint methods 
can be extended to more complex problems. 

In practice, flows are non-parallel and often three-dimensional. Compressibility 
may also be significant in some cases. How are receptivity characteristics to be 
found for such flows? Recent developments in the use of Parabolised Stability 
Equations (PSE) offer a promising possibility. By formulating and solving a set 
of adjoint parabolised equations, we have developed a method for mapping the 
efficiency with which external forcing excites the three-dimensional motions of a 
non-parallel boundary layer. The method makes use of the same computationally 
efficient formulation that makes the PSE currently so appealing. 

In the area of flow control, adjoint systems offer a powerful insight into the effect 
of control forces (Hill 1992). One of the simplest control strategies for boundary 
layers involves the application of localized mean wall suction. Why does it work 
so well? The adjoint method reveals a very simple flow analogy and a concise 
description of the effect of mean localized suction. 
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2. Accomplishments 

There are four areas where progress has been made. Firstly, the response of a 
two-dimensional incompressible parallel (Blasius) boundary layer to direct forcing 
has been investigated. This defines the elementary receptivity characteristics of a 
boundary layer. Secondly, a variety of natural forcing problems have been solved 
in which a scattering agent, such as surface roughness, couples freestream acoustic 
waves to Tollmien-Schlichting waves. Direct forcing of a non- parallel boundary layer 
is the third topic: here the adjoint to the Parabolised Stability Equations (PSE) 
is employed to deal in a computationally efficient manner with the non- parallel 
aspects of the problem. Finally, in the area of flow control, a new perspective is 
offered on the controlling effect of localized mean wall suction. 

2. 1 Direct forcing of the Blasius boundary layer 

In last year’s annual research brief, it was reported that the eigensolutions of 
the adjoint Orr-Sommerfeld equation, when suitably normalized, provide a detailed 
description of the response of a boundary layer to direct forcing. The characteristics 
of the adjoint to the Tollmien- Schlichting wave have been investigated, thereby 
developing a picture of the elementary processes by which Tollmien- Schlichting 
waves are produced most effectively. Software has been developed to determine the 
necessary normalized eigensolution of the adjoint Orr-Sommerfeld equation. 

The most significant features of the adjoint eigensolution, and consequently the 
physical properties of the boundary layer when subjected to direct time-harmonic 
forcing (i.e. an external source), are summarized a s follows: 

1. Over a wide range of frequencies and Reynolds numbers, the adjoint stream func- 
tion corresponding to the Tollmien-Schlichting eigensolution has a simple maximum, 
and far from the wall, it decays exponentially. 

2. The boundary layer is most sensitive to streamwise forcing (a momentum source) 
in the vicinity of the critical layer — the height above the wall at which the flow 
speed and the phase speed of the Tollmien-Schlichting wave coincide. The most 
sensitive y-location is shown in Figure 1 as a function of Reynolds number R — 
y/UooL/v , and frequency / = 27r/V/i7^ 0 , Here, L is the distance from the leading 
edge of the plate, is the flow speed at infinity, v is the viscosity, and /* is the 
frequency in Hertz. The solid line indicates the height above the wall at which 
streamwise forcing is most effective (the location of the maximum of the adjoint 
streamwise velocity component). The dashed line defines the position of the critical 
layer. 

3. Forcing in the wall-normal direction is much less effective than forcing in the 
streamwise direction. 

4. At the wall, normal motions create Tollmien-Schlichting waves much more effec- 
tively than do streamwise motions. 

5. The amplification of the Tollmien-Schlichting waves as they travel through the 
unstable region dictates that forcing at streamwise positions close to the lower 
branch leads to the strongest response. 
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FIGURE 1 . Height above wall at which a Blasius boundary layer is most sensitive to 
direct streamwise forcing, as a function of Reynolds number, for various frequencies 
/ x 10 6 . Dashed lines define the location of the corresponding critical layer. 


Detailed results are reported in Hill (1993). 

2.2 Natural forcing of the Blasius boundary layer 

Sound waves in the free stream can be scattered strongly into Tollmien- Schlicht- 
ing waves if there is even a weak mean flow distortion containing lengthscales com- 
mensurate with those of the Tollmien- Schlicht ing waves. The flow distortion might 
be caused typically by surface roughness or by mean suction at the wall. 

The ‘natural response’ problem in which a small surface roughness element acts as 
a scattering agent has been investigated by several researchers. Goldstein ( 1983) and 
Ruban (1985) used triple-deck theory to analyze the process in the infinite Reynolds 
number asymptotic limit. Crouch (1992) and Choudhari & Streett (1992) provide 
a solution of the incompressible problem at finite Reynolds number. The solutions 
indicate that the amplitude of the Tollmien-Schlichting wave that is produced by 
this mechanism is given by the product of an efficiency factor, a geometry factor, 
and the amplitude of the acoustic wave. The efficiency factor is a complex constant 
depending on the frequency and Reynolds number, and the geometry factor is the 
Fourier transform at the Tollmien-Schlichting wavelength of the roughness shape. 
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In the absence of any roughness, there is a profile U(y) on top of which is super- 
imposed an unsteady motion 

v(r,t) = u(y)e~ ,UJl x, where u(y) = 1 - e _(u ' fl/r2) ' (1) 

The planar fluctuations represent a freestream acoustic wave of unit amplitude 
which has a Stoke’s wave signature close to the plate. The frequency w is defined 
as f/R, and a is the TS waveumber at that frequency and Reynolds number. 

The roug hn ess modifies the mean flow, and the interaction of the Stoke’s wave 
with this mean flow distortion produces a TS-wave. It is assumed that, far from 
the roughness patch, the flow field recovers sufficiently quickly that the scattering 
takes place in an interaction zone in the vicinity of the roughness. 

The efficiency factor for the scattering process is found by examining the solution 
to the following inhomogeneous adjoint problem for the stream function 
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where ( u aul x + v au/ y)e~ ,iax ~ u,t) is the normalized adjoint eigensolution correspond- 
ing to the TS-wave. 

The solution 4 has some useful properties. The amplitude of the instability 
induced by the scattering of the freestream disturbance is 


A h(a) , 


where k(a) = 



(4) 


and the efficiency factor 


1 / dU du du a ^ \ . 

R\ dy dy 2 dy dy ) y=o 

The equations (2, 3) have been solved numerically, and A evaluated. The results 
are identical with those of Crouch (1992) and Cougar & Streett (1992). 

The solution to (2, 3) can also be used directly in the configuration in which mean 
suction at the wall acts as the scattering agent. Consider a velocity distribution 
V_ s (x) representing a suction/blowing distribution on the plate. The amplitude of 
the Tollmien-Schlichting wave is 


/ i d 3 4 . 

\iail dy 3 ^ 



V s (x)e~ ia *dx. 


( 6 ) 
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Surface admittance is defined as the ratio at the wall of the unsteady normal 
velocity to the unsteady pressure, and thus can be used to represent how the surface 
responds dynamically to unsteady pressures. Since there is no distortion of the 
mean flow, in this case the solution 4 is not employed. Suppose that there are 
spatially-uniform fluctuations p 0 e~ ,LJt in the pressure field. If 0 w (x) is the surface 
admittance, a Tollmien-Schlichting wave of amplitude 


/: 


|paw(0)| Po / Pw(x)e lQZ dx y 


( 7 ) 


will be induced. The adjoint pressure at the wall, p QW (0), associated with the 
normalized adjoint eigensolution in this case defines the efficiency factor for the 
scattering of freestream pressure fluctuations into Tollmien-Schlichting waves. 

2. S Direct forcing of non-parallel flow 

There has been considerable development in recent years in modeling transitional 
flows by the use of Parabolised Stability Equations (Herbert & Bertolotti 1987, 
Bertolotti 1991). The linear and non-linear dynamics of convectively-unstable dis- 
turbances in spatially-evolving boundary layers can be described accurately with 
little computational effort. The flow can be compressible, and the disturbances 
three dimensional. 

As the name suggests, the method involves “parabolising” the governing equa- 
tions for boundary layer disturbances. The solution is represented by a disturbance 
pattern resembling the local eigenfunction, modulated by a spatially-evolving oscil- 
latory factor. Both the disturbance pattern and the wavelength of the oscillation are 
assumed to evolve slowly with streainwise position. Starting at a chosen streamwise 
station, the solution is marched downstream; with a single sweep the evolution of 
the boundary layer disturbance is described. 

The formulation and solution of a set of adjoint Parabolised Stability Equations 
promises to provide a description of the efficiency with which a wide range of bound- 
ary layer motions are excited by direct forcing. In contrast with the regular PSE, 
the adjoint equations Eire marched upstream, starting at the outflow end of the com- 
putational domain. In this way, the events within the domaiin that give rise to a 
Tollmien-Schlichting disturbance at the outflow are identified. This approach is an 
extension to non-parallel flows of the results described in Section 2.1. The adjoint 
PSE can be solved within the same well-established computational framework as 
the regular PSE. The adjoint solutions are a natural complement to the regular 
solutions. 

Thanks to F. Bertolotti, a copy of the PSE library of subroutines has recently 
been made available. The following preliminary results have been obtained: 

1. The adjoint Parabolised Stability Equations have been formulated for three- 
dimensional disturbances in a two-dimensional spatially-evolving boundary layer. 

2. Solutions for a two-dimensional Tollmien-Schlichting wave in parallel flow have 
been checked. Figure 2 gives a graphic illustration of the receptivity maximum that 
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FIGURE 2. Plot of the efficiency with which Tollmien-Schlichting waves are excited 
by a point source of momentum of unit magnitude, oscillating at a frequency / = 
20 x 10 -6 . Contour values are shown in order to indicate relative magnitudes only. 

appears in the Blasius boundary layer (the y-scale has been expanded for purposes 
of visualization). It is positioned at the lower branch of the neutral stability curve 
and at a height of about half the displacement thickness from the wall, i.e. at 
the critical layer. (Figure 2 is a plot of the magnitude of the adjoint velocity 
and represents the magnitude of the response due to a unit amplitude harmonic 
mr.mpnt. nm source.) The solution has not been normalized, so that contour values 
do not indicate a physical measure of the response that will arise for unit forcing. 
However, the relative magnitudes illustrate the zone of high sensitivity. 

The effect of non-parallelism is expected to play a larger role for three-dimensional 
disturbances and compressible flows. This has yet to be investigated. 

2.4 Boundary layer control by auction 

Small amounts of localized wall suction can reduce dramatically the amplitude of 
Tollmien-Schlichting waves travelling in a boundary layer. This significant effect has 
been studied in detail because of its impact on Laminar Flow Control technology 
(Nayfeh at al. 1986, Saric k Reed 1986, Reynolds k Saric 1986). The numerical 
perturbation scheme of Reed k Nayfeh (1986) provides a computational analysis of 
the effect of an arbitrary distribution of suction strips beneath an incompressible 
boundary layer. Masad k Nayfeh (1992) have developed a scheme for compressible 
boundary layers. 

The effect of suction is to modify the mean flow both upstream and downstream 
of the slot. A TS-wave that enters this region of the flow grows at a rate which 
is different from that in the undisturbed flow. Integrating these changes over the 
entire flow gives the net effect of the suction on the disturbance amplitude. 

We construct here an analogous flow, i.e. a flow which has an identical effect 
upon the TS-wave amplitude. The analogous flow involves a local modification to 
the boundary layer profile directly above the slot, in proportion to the amount of 
suction at that streamwise station. 
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Figure 3. Schematic representation of the real flow (a) and analogous flow (b). 

The analogous flow replaces the boundary profile by a faster, thinner profile. A 
simple explicit form has been found for the equivalent profile: Let <p(x) be the 
strength of the wall suction velocity (scaled on Uoo) at streamwise location x , with 
local Reynolds number R. Let U(y) is the form of the profile, y being the distance 
from the wall. The net effect of the suction <f>(x) is identical to that of replacing 
U(y ) locally by 

Ua(y) = U(y) + 4>(x)R / ( 1 - U(y))dy, (8) 

Jo 

together with a uniform downflow of strength 4>(x) across the entire boundary layer. 

The fractional increase in the effective freestream flow speed is the product of 
<j>{x) with the local Reynolds number based on displacement thickness. 

Figures 3(a) and 3(b) offer a schematic representation of the real flow and the 
analogous flow, respectively. If the effect of the modification to the local growth 
of a TS-wave is integrated in each case, the same net change in amplitude will be 
found. 

To understand the effect of suction, we can thus consider the dynamics of Tollmien- 
Schlichting waves in thinner, faster boundary layers. For frequencies close to the 
lower branch of the neutral stability curve, the destabilizing influence of an increase 
in the flow speed is too weak to counter the stabilizing effect of the thinning of the 
layer. In practice, a disturbance at the frequency which is most “dangerous” from 
the point of view of transition is controlled by suction applied close to the lower 
branch. The disturbance amplitude is reduced typically by a large amount. By 
contrast, for frequencies close to the upper branch (much less “dangerous” from the 
point of view of transition), the modified profile tends to be less stable. Suction 
leads to an increase in the amplitude of disturbances at these frequencies. However, 
there is no reported experimental evidence of amplification of higher frequencies 
due to suction. 

3. Future work 

The efficiency with which acoustic waves are scattered into TS-waves by surface 
roughness hits already been investigated in detail. The acoustic waves are taken 
to have infinite wavelength, which reflects the “infinitely” fast speed at which they 
propagate. Disturbances such as freestream turbulence convect at, or close to, the 
flow speed and have finite lengthscales associated with them. How efficiently do 
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such motions scatter into TS-waves? I intend to address this question using an 
extension of the adjoint method employed for the acoustic scattering problem. 

Secondary instabilities play a key role in the so-called K- and H-type transition 
routes: three-dimensional disturbances grow upon a finite amplitude TS-wave until 
the flow evolves rapidly to turbulence. The following questions will be addressed: 
What is the most efficient means of exciting the secondary instability? Can it be 
controlled /suppressed? 

Stationary crossflow vortices appear on swept airfoils in response to surface rough- 
ness. Secondary instabilities then lead to a breakdown to turbulence. It is impor- 
tant to understand the process by which the crossflow vortices arise and to identify 
those locations where roughness elements are most important. An investigation of 
this problem will be made using the adjoint PSE in combination with the classical 
independence principle. 
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